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ABSTRACT 

During their infancy, stars are well known to expel matter violently in the form of well-defined, 
collimated outfiows. A fairly unique exception is found in the Orion BN/KL star- forming region 
where a poorly collimated and somewhat disordered outfiow composed of numerous elongated "finger- 
like" structures was discovered more than 30 years ago. In this letter, we report the discovery in the 
same region of an even more atypical outfiow phenomenon. Using CO (2-1) line observations made 
with the Submillimeter Array (SMA), we have identified there a 500 to 1,000 years old, expanding, 
roughly spherically symmetric bubble whose characteristics are entirely different from those of known 
outflows associated with young stellar objects. The center of the bubble coincides with the initial 
position of a now defunct massive multiple stellar system suspected to have disintegrated 500 years 
ago, and with the center of symmetry of the system of molecular fingers surrounding the Kleinmann- 
Low nebula. We hypothesize that the bubble is made up of gas and dust that used to be part of 
the circumstellar material associated with the decayed multiple system. The Orion hot core, recently 
proposed to be the result of the impact of a shock wave onto a massive dense core, is located toward 
the south-east quadrant of the bubble. The supersonic expansion of the bubble, and/or the impact of 
some low-velocity filaments provide a natural explanation for its origin. 

Subject headings: stars: formation - ISM: jets and outfiows - techniques: imaging spectroscopy - 
techniques: high angular resolution - ISM: molecules 



1. INTRODUCTION 

Discovered more than 40 years ago, the Becklin- 
Neugebauer /Kleinmann-Low ( BN/KL) region in Orion 
(jBecklin fc Neugebauerl 119671 : iKleinmamT Lowl [19671 ) 
remains quite enigmatic. It is at the center of a 
fast (30-100 km s~^), and massive (^ 10 M©) out- 
fiow particularly prominent in mole cular hydrogen emis- 
sion at near-infrared wavelengths (Kwa n fc Scoviiyil976l : 
I Allen fc Burton 1993: Tavlor et al. 1981 ). Composed of 
numerous elongated "finger-like" structures, this out- 
fiow has a very peculiar morphology suggestive of an 
explosive origin. A strong and turbulent hot molecu- 
lar core (Wright et al.l ll996l : iBlake et al.lll996l : iLiu et al.l 
[2002) associated with i ntense maser emiss i on from sev- 
eral mole cular spe cies (|Gaume et al.l 1 19981 : iCohen et al.l 
120061 : [Zapata et aDl2009aD also lies in the BN/KL region, 
but the very nature of this hot core and its relation to 
the outfiow remain debated. Three massive young stars 
(BN, / and n) are known to be moving away from the KL 
nebula (Plambeck et al.l 119951 : IGomez et al.ll2005L I2008D 
at a few tens of km s~^ . Their velocity vectors point away 
from a common point of origin where the three stars must 
have been located about 500 years ago. This suggests 
that they were originally members of a multip le system 
that dynamica lly decayed very recently (Rodriguez et al. 
120051 : IGomez et al...2005. .2008 ). Remarkably, our recent 
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molecula r observations of th e outflow associated with 
BN/KL (jZapata et al.ll2009b[ ) showed that the elongated 
fingers forming that flow also point exactly away from 
the position where the initial multiple stellar system was 
located. Moreover, in this study, we found that the dy- 
namical age of the outflow is consistent with about 500 
years. It is, therefore, very likely that the same energetic 
phenomenon that led to the acceleration of the stellar 
sources BN, / and n was also responsible for the explo- 
sion that gave birth to the peculiar gaseous outfiow in 
the BN/KL region of Orion. 

2. OBSERVATIONS 

In order to further explore the interrelation between 
the various energetic phenomena occurring in Orion 
BN/KL, we have carried out interferometric observations 
of the region in the optically thin 2^1 line of ^^CO at 
u = 220.3 GHz using the Submillimeter Array (SMA) 
()Ho et al.ll2004f ). This transition is a reliable indicator 
of the kinematics and mass of the region observed. The 
observations were collected in 2007 January and 2009 
February while the SMA was in its compact and sub- 
compact configurations, with baselines ranging in pro- 
jected length from 6 to 58 m; this provides an angular 
resolution of about 4". The digital correlator was config- 
ured to provide a spectral resolution of 0.40 MHz, cor- 
responding to a resolution in velocity of 1.05 km s~^. 
The zenith opacity {r230GHz) was ~ 0.1-0.3, indicating 
reasonable weather conditions. Observations of Uranus 
and Titan provided the absolute scale for the flux den- 
sity calibration. The gain calibrators were the quasars 
0530+135, 0541-056, and 0607-085. 

The data were calibrated using the IDL superset MIR, 
originally developed for the Owens Valley Radio Obser- 
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Fig. 1. — Color and contour images of the blueshifted (upper panel) 
and redshifted (lower panel) integrated "'^'^00(2-1) line emission from 
the molecular bubble overlaid with the ""^^00(2-1) filaments. Upper 
panel: The white contours represent the blueshifted ^'^CO(2-l) emis- 
sion and are at —2, 2, 3, and 4 times 2.5 Jy beam~^ km s"""^. The blue 
contours represent blueshifted ^^CO(2-l) emission and are at —4, 4, 8, 
16, 24, 32, 38, 42, and 52 times 2 Jy beam~^ km s~^ Lower panel: The 
yellow contours represent the redshifted ^"^CO(2-l) emission and are at 
— 2, 2, 3, 4 and 5 times 13 Jy beam~^ km s~^. The orange contours 
represent redshifted ^^CO(2-l) emission and are at —4, 4, 8, 16, 24, 32, 
38, 42, and 52 times 2.7 Jy beam~^ km s~^. In both panels the blue 
and red points represent the position of the ""^^00(2-1) emission peaks 
in every velocity channel. The synthesized beam of the line image is 
4'.'3 X 3'.'2 with a P. A. of -4.0° and is shown near the lower left corner of 
each panel. The wedge indicates the line emission in mJy beam~^ km 
s~^. The red and gr een filled dots mar k the positions of the runaway 
sources 57V, / and n (jGomez et al.|[2'QQ5r) . The v ectors on these sources 
represent the direction o f their proper motion (iRodriguez et ani20Q5l ; 
iGomez et al .''2005. 2008"). The red and green cross-hair circles marks 
the zone from where the three sources were ejected some 500 years ago 
(IGomez et al.ii200 5) and the origin of the explosive flow (Zapata et al.' 
|2009bD . The dashed lines in the top panel mark the position and ori- 
entation of the position-velocity diagrams presented in Figure 2 were 
computed. 



vatory CScovi lle et allll993 r) and adapted for the SMAEI- 
The cahbrated data were imaged and analyzed in a stan- 
dard manner using the MIRIAD, GILD AS, and AIPS 
packages. We used the ROBUST parameter set to 2 to 
obtain a slightly better sensitivity sacrificing some angu- 
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Fig. 2. — ^"^CO(2-l) position-velocity digrams of the molecular ring 
computed at a position angle of +210° (upper panel) and +300° (lower 
panel). The contours in both panels are from 26%, 27%, 29%, 30%, 
40%, 50%, 60%, 70%, 80% and 90% of the peak of the line emission. 
The systemic LSR radial velocity of the ambient molecular cloud is 
about 9 km s~^. The synthesized beam of the line image is 4'.' 3 X 3'.' 2 
with a P. A. of —4.0°. The pink star marks the position and LSR ve- 
locity of the center of the explosive outflow (Zapata et al. 2009b). The 
emission at ambient velocities is clearly extended and poorly sampled 
with the SMA. In the right bottom corner of the upper panel, we show 
a intensity-velocity cut of the ring at a P. A. = +90° (see the white 
line) which reveals that expansion velocity of bubble is approximately 
15 km s"""^. In the bottom panel, we mark the position of the Orion 
Hot Core. 

lar resolution. The line image rms noise was around 100 
mJy beam~^ for each channel. 

3. RESULTS AND DICUSSION 

The CO (2-1) line detected in our observations is 
broad, extending from about —30 km s~^ to +40 km s~^ 
in radial velocity. The structure of the emitting region 
is shown in Figure 1, where we present separately the 
blueshifted (from —31 to —18 km s~^; upper panel) and 
redshifted (+23 to +33 km s~^; lower panel) integrated 
intensity maps. In Figure 1, we also show (as contours) 
the integrated intensity of the high- velocity fingers (from 
— 130 to —35 km s~^ for the approaching velocities, and 
from +35 to +120 km s~^ for the receding velocities) 
traced by ^^C0(2-l) emission (jZapata et al.ll2009bf ). and 
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the positions and proper motions of the runaway sources 
(BN, / and n). The position of the multiple stellar sys- 
tem which contained these sources before it decayed is 
shown as a crossed-haired circle, and we shall henceforth 
refer to this point as "O" . 

The CO (2-1) emission clearly traces a roundish 
structure curling around position "O". The emission at 
position "O" itself is weak so the structure must be hol- 
low: it is either a ring or a shell. In Figure 2, we show 
the CO (2-1) position- velocity diagrams of that struc- 
ture computed along strips centered on "O" along posi- 
tion angles of +210° and +300°. The morphology seen 
in these diagrams is characteristic of a bubble: redshifted 
and blueshifted emission are simultaneously present to- 
ward the center of the structure, corresponding to the 
two sides of the bubble. It is impossible from its mor- 
phology and kinematics alone to decide whether the bub- 
ble is expanding or contracting. The fact that the three 
stars (BN, /, and n) known to be moving away from 
point "O" are projected exactly on top of the bubble, 
however, almost certainly indicates that it is expanding. 
The expansion velocity of the bubble is of order 15 km 
(Fig ure 2) and, given the distance to Orion (414 pc; 
iMenten et al.l 120071 ). its radius is about 2000 AU. This 
yields a dynamical age for the bubble of about 600 years, 
with an uncertainty of a few hundred years. 

The C^^ 0(2-1) line was also included in our SMA ob- 
servations. The bubble is marginally seen in these data, 
with a typical intensity 6 to 8 times weaker than in the 
i3CO(2-l) line. This indicates that the ^^00(2-1) line 
is optically thin, and can be used to estimate the mass 
of the expanding bubble. Assuming an excitation tem- 
perature of 50 K, and an abundance ratio between II2 
and ^^CO of 5 x 10^, we obtain a mass of about 5 Mq 
for the entire bubble. The corresponding kinetic energy 
(assuming an expansion velocity of 15 km s~^) is about 
1 X 10^^ erg, roughly one order of magnitude smaller 
than the kinetic energy of the outflow (iKwan fc Scovilld 
[Hzi). 

The bubble is centered on point "O" , which decidedly 
appears to be a very special position within the Orion 
BN/KL region. This, combined with the remarkably 
similar timescales of all three phenomena, clearly indi- 
cates that the acceleration of the massive young stars 
BN, /, and n, and the births of both the fast outflow and 
the bubble most certainly occurred as a consequence of 
a single energetic event. Indeed, a natural scenario ex- 
plaining the creation of the bubble follows from the in- 
terpretation of the large velocities of BN, /, and n in 
terms of a dynamical disintegration. Young stars are 
well-known to be surrounded by circumstellar disks and 
envelopes. If a multiple system disintegrates because of 
a chaotic few-body encounter, the circumstellar mate- 
rial initially in those structures will clearly be affected. 
The gas most strongly bound to the young stars may be 
ejected together with the stars themselves, and some of 
the material may even be independently ejected to high 
speed by a sling-shot effect. The least strongly bound 
material, however, will likely remain nearly stationary 
at first. However, since the stars that provided most of 
the gravitational binding energy of the young system are 
now gone, this material will find itself with an excess of 
kinetic energy, and will start to expand at a speed typical 



of the velocity dispersion of the material (^15 km s~^ as 
we will see momentarily). Moreover, since it was initially 
distributed among several disk/envelope structures with 
randomly oriented angular velocities, this material would 
be expected to expand fairly isotropically rather that as 
a highly flattened system. Finally, because of conserva- 
tion of angular momentum, the system will progressively 
loose any original flattening it might have had. 

As mentioned earlier, the kinetic energy of the bubble 
is about 1 X 10^^ erg, a factor of 100 to 1000 times less 
than the disruptive energy of /, and n and BN 500 years 
ago (|Ballv fc Zinneckerll2005l ). Its gravitational energy, 

on the other hand, is about = — 2 x 10^^ erg. 

Thus, the bubble is clearly not currently in virial equi- 
librium, and must be expanding -in agreement with our 
previous considerations. Before the disintegration, how- 
ever, BN, /, and n as well as the 15 Mq in the outflow 
and bubble were all part of a common system. The com- 
bined mass of BN, /, and n must be of order 35 M©, 
so the total mass of the system would have been about 
50 Mq. For such a system to be in virial equilibrium, 
it must be confined to a radius of about 200 AU. Inter- 
estingly, this is approximately the size of the region (1^^) 
where the sources BN, /, and n were located about 500 
years ago (iGomez et al.l [2QQ5I ) . Moreover, the typical 
speed of the gas in such a virialized structure would be 

V = \J~^^ ~ 15 km s~^, in excellent agreement with the 

measured expansion velocity of the bubble. 

The expansion velocity of the bubble is clearly super- 
sonic, so shocks should be created where it encounters 
dense gaseous clumps during its expansion. Interest- 
ingly, our SMA observations do reveal the presence of 
typical shock tracers {e.g. SiO, SO, and SO2) in the 
bubble. The Orion hot core mentioned earlier is located 
toward the south-east quadrant of the bubble (see Fig- 
ure 2), and was recently suggested to be caused by the 
impact of a shock wave related to the dynamical decay 
of BN, /, and n onto a dense pre-existing core called the 
"Extended Ridge" (Zapata et al. 2010). To explain the 
near systemic mean velocity of the hot core, the authors 
of that suggestion argued that the shock wave must be 
fairly slow, and i nteract with a part icularly dense and 
massive material. (jZapata et al.|[2QTQ ^ proposed that the 
required shock wave might be due to the impact of low- 
velocity filaments. However, the supersonic expansion 
of the bubble reported here might provide a more nat- 
ural explanation of the shocks required to explain the 
Orion Hot Core. For a more complete discussion of the 
nature of the Orion Hot Core, we refer the readers to 
IZapata et aH (j2010[ ). 

4. CONCLUSIONS 

In this letter, we reported the discovery of an entirely 
new type of outflow phenomenon in the BN/KL region of 
Orion. Using CO (2-1) SMA observations, we identified 
a 500 to 1000 years old, roughly spherically symmetric 
bubble, expanding at about 15 km s~^. The center of 
the bubble coincides with the initial position of a now de- 
funct massive multiple stellar system suspected to have 
disintegrated 500 years ago, and with the center of sym- 
metry of the system of molecular fingers surrounding the 
Kleinmann-Low nebula. We propose that the material 
in that bubble was originally associated with the young 
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stars in the now defunct multiple system. 
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